Amorphization induced by pressure: results for zeolites and general implications 
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We report an ab initio study of pressure-induced amorphization (PIA) in zeolites, which are model 
systems for this phenomenon. We confirm the occurrence of low-density amorphous phases like the 
one reported by Greaves et al. [Science 308, 1299 (2005)], which preserves the crystalline topology 
and might constitute a new type of glass. The role of the zeolite composition regarding PIA is 
explained. Our results support the correctness of existing models for the basic PIA mechanim, but 
suggest that energetic, rather than kinetic, factors determine the irreversibility of the transition. 

PACS numbers: 61.50Ks,61.43.Fs,64.70.Pf,82.75.-z 
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Many minerals can be turned amorphous by mere ap- 
plication of pressure, an intriguing phenomenon known 
as pressure induced amorphization (PIA) Q]. The con- 
crete realizations of PIA can be quite diverse. Remark- 
ably, upon release of the applied pressure some materials 
recover their crystalline order, thus exhibiting reversible 
amorphization, while others remain amorphous. Our un- 
derstanding of PIA is only partial. Simulations of sys- 
tems like a-quartz 0, silica and a-berhnite 0| sug- 
gest PIA is the result of a first-order transition associated 
with very localized, weakly interacting structural distor- 
tions that become unstable upon compression. In such a 
transition, domain nucleation would overwhelm growth 
and destroy the long-range order . Cohen, Iniguez and 
Neaton [| (CIN) further propose the PIA transition will 
be reversible if the crystalline topology is preserved in 
the amorphous phase, i.e., if the amorphization does not 
involve bond formation or breaking. While physically 
plausible, this picture is yet to be confirmed. 

The recent work of Greaves et al. on the nanoporous 
aluminosilicates known as zeolites has renewed the in- 
terest in PIA. These authors have shown a zeolite 
may present two distinct PIA phases: a low-pressure 
(« 2 GPa) low-density amorphous phase (LDA), which 
they argue may constitute a new type of glass, and a 
high-pressure (« 6 GPa) high-density amorphous phase 
(HDA). Further, they claim the crystalline topology is 
preserved in the LDA phase and lost in the HDA phase, 
which, according to the CIN picture, implies amorphiza- 
tion will be reversible in the former case and irreversible 
in the latter. These results, together with other studies 
0,0 that, for example, show a striking dependence of the 
PIA reversibility on the zeolite composition, clearly point 
at these systems as ideal to test general PIA theories. 

Here we present an ab initio study of PIA in three 
representative zeolites with different compositions. Our 
results (i) confirm the occurrence of the above mentioned 
LDA and HDA phases, (ii) show how the zeolite compo- 
sition controls the nature of the PIA transition, and (iii) 
essentially confirm, and perfect, the CIN picture regard- 
ing PIA reversibility. While obtained for zeolites, our 



results clearly pertain PIA phenomena at large. 

We used the Generalized Gradient Approximation to 
Density Functional Theory as implemented in the 
code SIESTA |0. Note that we wanted to study amor- 
phization occurring in spite of neglegible thermal ac- 
tivation, and thus focused on low temperature simula- 
tions. We proceeded as follows: at each considered pres- 
sure, we performed a short (100 fs) molecular dynam- 
ics simulation at 100 K, starting with random velocities, 
and relaxed the resulting structure. Zeolites have the 
general formula ^"^^Ala;Sii_a;02, where ^ is a charge- 
compensating cation (e.g. Li or Na). The Si and Al 
atoms are at the center of corner-sharing O4 tetrahedra. 
A typical zeolite structure, with the so-called LTA frame- 
work, is sketched in the top-right inset of Fig. ^ note 
the four-, double- four-, six-, and eight-member rings (de- 
noted as 4MR's, D4R's, etc.), also known as secondary 
building units. We studied three LTA zeolites with differ- 
ent percentages of Al and Na as the charge-compensating 
cation: an Al-free "all-Si02" system, Na-ZK4 with a 1-to- 
5 Al-Si ratio, and Na-A where the ratio is 1-to-l. Follow- 
ing experimental information [ll] |. we considered primi- 
tive cells containing 72, 76, and 168 atoms, respectively, 
for all-Si02, Na-ZK4, and Na-A. 

All-Si02 results.- The all-SiOa LTA zeolite exhibits 
most of our key findings. Figure ^ shows the pressure 
dependence of the unit cell volume and the evolution of 
an 8MR that captures the typical structural distortions. 
The slope changes and volume discontinuities indicate a 
series of phase transitions. It is sufficient for our pur- 
poses to describe the structure in terms of the average 
values and standard deviations of the relevant angles (Si- 
0-Si and 0-Si-O) and distances (Si-0). For example, 
for the reference LTA structure at GPa we obtained 
Jsio=l-63±0.01 A, 0osiO=lO9±l°, and ^siosi=151±8°. 

Phases I to III in Fig. ^ are connected by continu- 
ous transitions characterized by rigid rotations of the 
O4 tetrahedra. Such Rigid Unit Modes (RUM's 0) 
are mainly reflected in ^siOSi- For example, for 
phase III at 2.75 GPa we obtained dsio=1.62±0.01 A, 
^osiO=109±2°, and 6isiosi=142±18°. At 3.5 GPa 
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phase III transforms discontinuously into a phase IV of 
significantly smaUer volume. This transition does not in- 
volve any topological change, i.e. no bonds are formed 
or broken. In fact, as in the previous cases, the struc- 
tural changes mostly affect ^siosi; at 7 GPa we obtained: 
dsio=l-63±0.02 A, ^osio=109±7°, and 0siOSi=132±2O°. 

At 11.25 GPa phase IV transforms discontinuously into 
a phase V, and at 12.25 GPa there is another first- 
order transition to a phase VI. Phases V and VI dis- 
play collapsed rings of all types and new Si-0 bonds 
(see the 8MR depicted in Fig. with the corre- 
sponding loss of the LTA-framework topology. The oc- 
currence of SiOs and SiOg groups results in a wide 
dispersion of distances and angles. For example, at 
13 GPa we obtained Jsio=l-70±0.10 A, ^osiO=107±23°, 
and 6'siosi=117±17°. Note the increase in the aver- 
age Si~0 distance, which reflects the existence of high- 
coordination defects. 

FigureHalso shows our results regarding the reversibil- 
ity of the transtions. The transitions to phases II and 
III are reversible, with no hysteresis in the V{p) curve. 
Phase IV can also transform back to the crystalline 
phase, but hysteresis occurs in this case. Note that the 
presence or absence of hysteresis agrees with the observed 
character, first- or sencond-order, of the transition. Fi- 
nally, upon decompression from phases V and VI, the sys- 
tem undergoes a number of structural changes but does 
not find its way back to the low-pressure stable phases. 
The transitions to phases V and VI are thus irreversible. 

Our simulations reveal the atomistic origin of this irre- 
versibility. The metastable phases V** and VP* in Fig.Q] 
present the right first-neighbor coordination and can thus 
be viewed as formed by Si04 units. However, the ring 
structure is not the LTA one. For example, in phase V** 



FIG. 2: All-Si02 2x2x1 supercells (see text) resulting from 
transitions at 3.50 GPa (panel a) and 11.25 GPa (panel b). 
Dashed lines sketch the cell we start from. We highlight ex- 
amples of atoms that were translationally related before the 
transition; a SiOs group is circled in panel b. 

the original 8MR has transformed into two disconnected 
4MR's, and phases V* and VI** contain Si pairs shar- 
ing two oxygens (see Fig.^. Note that phases V** and 
VI** are robustly metastable: while their excess energy 
with respect to phase I is relatively large (about 0.35 eV 
per formula unit at GPa) , a transition to the crystalline 
phase would require multiple Si-0 bond breakings within 
Si04 units, which is energetically very costly. 

The origin of the defects affecting the ring topology can 
be easily identified. Phases V and VI exhibit collapsed 
rings in which bonds form between Si and O atoms on op- 
posite sides of the original rings. The resulting SiOs and 
SiOe groups break upon decompression, and the atoms 
recover their original coordination. However, this defect 
breaking can happen in a variety of ways from which only 
one allows the ring to fully recover at low pressures. In 
addition, this unique right way involves the largest vol- 
ume expansion and, thus, is energetically favorable only 
at relatively low pressures. For example, for phase VI 
we obtained that above 2 GPa the SiOs and SiOe groups 
find it energetically favorable to break in ways that re- 
store the original atomic coordination but, at the same 
time, destroy the 8MR. 

We also investigated if these transitions involve true 
amorphization. Our simulated system is defined by the 
unit cell of the LTA structure and, thus, cannot cap- 
ture the translational symmetry breaking characterizing 
amorphization. Hence, in a few selected cases we consid- 
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FIG. 3: Same as Fig.[Tlfor the Na-A zeolite. 

ered larger supercells and simulated the loss of transla- 
tional symmetry directly. More precisely, we considered 
the relaxed structures obtained at pressures neighboring 
transition points (i.e. GPa, 3.5 GPa, etc.), created 
the corresponding 2 x 2 x 1 x 72-atom supercells, and in- 
creased the pressure to observe how the transformation 
proceeds. Interestingly, for the second-order transitions 
to phases II and III, the translations within the supercell 
are preserved, indicating that no PIA occurs. On the 
other hand, as shown in Fig. [21 the translational symme- 
try within the supercell is completely lost in the transi- 
tions to phases IV and V, both of which are first-order 
in character. We can thus conclude that phases IV and 
V are genuine amorphous phases, the crystal topology 
being preserved in the former and lost in the latter. 

Na-ZK4 and Na-A results - We consider first the case 
of Na-A, which is more informative. Figure 13 shows the 
pressure dependence of the 168-atom primitive cell and 
the evolution of a D4R, which captures the typical struc- 
tural distortions. (Note that the Si and Al atoms are 
intercalated in Na-A.) Phase I is the reference phase, 
which at GPa is characterized by JsiO = l-65±0.01 A, 
Jaio1-76±0.02 a, ^osio = 109±5°, ^oaio = 109±3°, and 
fi'siOAi=148±10°. Phases II and III are the product of 
two transitions that are, respectively, continuous and 
discontinuous. Both transformations are characterized 
by RUM's that involve rotations of the O4 tetrahedra 
and mainly affect the Si-O-Al angles. For example, 
for phase III at 5 GPa we obtained dsiO=l-66±0.05 A, 
dAio=l-82±0.09 A, 6'osiO=109±10°, 6Ioaio=108±21°, 
and 6'siOAi=118±18°. Note that phase III preserves 
the LTA topology, clearly resembling what we found for 
phase IV of all-SiOz. 

At about 5.25 GPa the system undergoes a second dis- 
continuous transition to a phase IV in which the D4R's 
collapse and the LTA topology is lost. Comparison with 
all-Si02 clearly suggests the presence of Al reduces the 
pressure at which the coordination defects appear, which 



is consistent with the fact that Al is more likely than Si 
to have a 5- or 6-fold oxygen coordination. Note also that 
only the D4R's collapse in phase IV. Interestingly, in Na- 
A the centers of all the 6 and SMR's are occupied by Na 
atoms, which suggests that the cations are responsible 
for the preservation of the large rings. The structure of 
phase IV at 7 GPa is characterized by Jsio=1.67±0.06 A, 
Jaio1-85±0.10 A, ^osiO=109±12°, ^oaio = 107±24°, and 
^SiOAi=116±20°. The relatively large standard devia- 
tions reflect the structural disorder. 

As shown in Fig. 13 we found that all the transitions in 
Na-A are reversible. Hysteresis does not occur in the case 
of the second-order transition (to phase II), but it does 
for the first-order transitions (to phases III and IV). Note 
that the reversibility from phases II and III, in which the 
ideal LTA topology is preserved, is consistent with our 
results for all-Si02. However, the obtained reversibil- 
ity from phase IV, in which the LTA topology is lost, is 
clearly at odds with what we found in the all-Si02 case. 
Note also that, unlike phases V** and VI** of all-Si02, 
phases IV* and IV** of Na-A present coordination defects 
and are not robustly metastable at low pressures. 

Our results for Na-ZK4 (not shown here) can be sum- 
marized as follows: (1) At low pressures there are contin- 
uous transitions dominated by RUM's. (2) A first-order 
transition, at 4.25 GPa, causes coordination defects, most 
of which involve neighboring Si-Al atom pairs that ap- 
proach to share two O atoms. All the Al atoms in the sys- 
tem give raise to coordination defects, while only a small 
fraction of the Si atoms do. (3) Upon further compression 
we find discontinuous transitions to phases in which rings 
of all types collapse. (4) The second-order transitions are 
reversible without hysteresis, and all the first-order tran- 
sitions are irreversible. Thus, as in Na-A, the presence 
of Al in Na-ZK4 reduces the pressure at which the coor- 
dination defects appear. On the other hand, at variance 
with all-Si02 and Na-A, Na-ZK4 does not present any 
first-order transition to a topology-preserving phase. 

Discussion.- Our simulations render a wealth of con- 
clusions that pertain not only zeolites, but PIA phe- 
nomena at large. Maybe most importantly, we confirm 
the existence of the LDA phases reported by Greaves 
et al. 0. Two of the considered zeolites (all-Si02 and 
Na-A) present an LDA phase, which supports the claim 
that such phases may be quite common. Further, for 
all-Si02 the LDA phase is predicted to be stable in a 
wide range of pressures. Our results thus suggest that 
this zeolite, which has been recently synthesized 0| and 
is relatively simple, would be ideal for detailed experi- 
mental and theoretical studies of the new type of glass 
proposed in Ref. 0. 

Our simulations confirm PIA is a first-order transi- 
tion, thus supporting the above mentioned nucleation- 
overwhelms-growth mechanism |^ . In topology-breaking 
PIA transitions, the structural distortions involve forma- 
tion of new bonds and, as Fig. [2b suggests, are rather 
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localized and can freeze in independently from the rest. 
In the PIA transitons that respect topology, the localized 
distortions are essentially rigid rotations of the O4 tetra- 
hedra. This is further confirmed by vibrational calcula- 
tions showing the whole low-energy RUM-related band 
softens under compression. 

The Al atoms facilitate the formation of coordination 
defects. Such defects are predicted to occur at about 
5 GPa in Na-ZK4 and Na-A, which is consistent with 
experimental data 001, and only above 11 GPa in all- 
Si02. The Na cations, on the other hand, impede the 
collapse of the rings at whose centers they are located. 
As a result, there is no formation of new Si-0 bonds be- 
tween atoms on opposite sides of Na-hosting rings. This 
seems consistent with our finding that Na-A, which con- 
tains a large amount of Na, presents an LDA phase while 
Na-ZK4 does not. Also, based on our results, it seems 
reasonable to assume large cations will be more effective 
in preventing rings from collapsing. That is in agreement 
with reports that, for one particular zeolite, PIA is irre- 
versible for small cations (e.g. H) and reversible for larger 
ones (e.g. Li and Na). The same rationale applies to the 
probable role of H2O molecules preventing PIA 

Our results support the GIN picture ^ that topology- 
preserving PIA transitions are reversible. However, at 
variance with what is proposed in Ref. 13, they also 
show that topology-breaking PIA transitions may be re- 
versible. That is the case of Na-A, where the topology 
breaking is caused by bond formation between atoms 
that are close neighbors in the crystalline phase. We 
find in such conditions the coordination defects can be 
correctly undone upon decompression, so that the crys- 
talline structure is recovered. We find the PIA transition 
to be irreversible in cases in which the coordination de- 
fects involve atoms that are away from each other in the 
crystalline phase (i.e. when the 6M and 8M rings col- 
lapse). In such cases, the atoms recover their preferred 
low-pressure coordination upon decompression, but the 
resulting ring topology differs from the crystalline one. 

Our simulations also refine the GIN picutre in what 
regards the mechanism for irreversibility. According to 
Ref. 0, irreversibility occurs because the system is un- 
able to find the transition path to its most stable (crys- 
talline) phase. To some extend our results support such 
a kinetics-related explanation, as we find that, for amor- 
phous phases with collapsed rings, only a small number 
of transition paths allow the system to recover both the 
ideal coordination to first neighbors and the LTA ring 
structure. However, our work also suggests a second and 
more important cause for the irreversibility: We find that 
at high pressures (e.g. above 2 GPa) it is energetically 
favorable to recover the ideal first-neighbor coordination 
in ways that break the crystalline ring structure. Thus, 
in addition to being more numerous, the wrong transi- 



tion paths (i) are energetically preferred at moderately 
high pressures and (ii) lead to phases that are robustly 
metastable at low pressures. This thermodynamics- 
related mechanism constitutes an alternative cause, prob- 
ably the main one in zeolites, for PIA irreversibility. 
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